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We study the interlayer coherent incompressible phase in Trilayer Quantum Hall systems (TLQH) 
at total filling factor vr = 1 from three approaches: Mutual Composite Fermion (MCF), Composite 
Boson (CB) and wavefunction approach. Just like in Bilayer Quantum Hall system, CB approach 
is superior than MCF approach in studying TLQH with broken symmetry. The Hall and Hall drag 
resistivities are found to be quantized at h/e 2 . Two neutral gapless modes with linear dispersion 
relations are identified and the ratio of the two velocities is close to v3- The novel excitation 
spectra are classified into two classes: Charge neutral bosonic 2-body bound states and Charge 
±1 fermionic 3-body bound states. In general, there are two 2-body Kosterlize-Thouless (KT) 
transition temperatures and one 3-body KT transition. The Charge ±1 3-body fermionic bound 
states may be the main dissipation source of transport measurements. The broken symmetry in 
terms of SU(3) algebra is studied. The structure of excitons and their flowing patterns are given. 
The coupling between the two Goldstone modes may lead to the broadening in the zero-bias peak 
in the interlayer correlated tunnelings of the TLQH. Several interesting features unique to TLQH 
are outlined. Limitations of the CB approach are also pointed out. 
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I. INTRODUCTION 

Quantum Hall Effect (FQHE) in multicomponent sys- 
tems has attracted a lot of attention since the seminal 
work by Halperin [1]. These components could be the 
spins of electrons when the Zeeman coupling is very small 
or layer indices in multi-layered system. In particular, 
spin-polarized Bilayer Quantum Hall (BLQH) systems 
at total filling factor vt = 1 have been under enormous 
experimental [2-5] and theoretical [6-13] investigations 
over the last decade. When the interlayer separation d is 
sufficiently large, the bilayer system decouples into two 
separate compressible v — 1/2 layers. However, when d is 
sufficiently small, in the absence of interlayer tunneling, 
the system undergoes a quantum phase transition into 
a novel spontaneous interlayer coherent incompressible 
phase [11]. 

By treating the two layer indices as two pseudo- 
spin indices, Girvin, Macdonald and collaborators [9-11] 
mapped the bilayer system into a Easy Plane Quantum 
Ferromagnet (EPQFM). They achieved the mapping by 
projecting the Hamiltonian of the BLQH onto the Lowest 
Landau Level (LLL) and then using subsequent Hartree- 
Fock (HF) approximation and gradient expansion ( called 
LLL+HF in the following ). They explored many rich 
and interesting physical phenomena in this system. The 
low energy excitations above the ground state is given 
by an effective 2 + 1 dimensional XY model. In addition 
to the Goldstone mode, there are also 4 flavors of topo- 
logical defects called " merons " which carry fractional 
charges ±1/2 and also have ± vorticity. They have log- 
arithmic divergent self energies and are bound into pairs 
at low temperature. The lowest energy excitations carry 
charge ±e which are a meron pair with opposite vor- 
ticity and the same charge. There is a finite tempera- 



ture Kosterlize-Thouless (KT) phase transition at Tkt 
where bound states of the 4 flavors of merons are broken 
into free merons. Unfortunately, this transition has not 
been observed so far. The EPQFM approach is a micro- 
scopic one which takes care of LLL projection from the 
very beginning. However, the charge sector was explicitly 
projected out, the connection and coupling between the 
charge sector which displays Fractional Quantum Hall ef- 
fect and the spin sector which displays interlayer phase 
coherence was not obvious in the EPQFM approach. 

In [12], I used and extended both Mutual Composite 
Fermion (MCF) and Composite Boson (CB) approaches 
to study both balanced and imbalanced BLQH and made 
critical comparisons between the two approaches. I iden- 
tified several serious problems with the MCF approach 
and then developed a simple CB approach which natu- 
rally avoids all these problems suffered in the MCF ap- 
proach. I found the CB approach is superior to MCF 
approach in the BLQH with broken symmetry. The func- 
tional form of the spin sector of the CB theory is the same 
as the EPQFM. It also has the advantage to treat charged 
excitations and the collective order parameter fluctua- 
tions in the pseudo-spin sector on the same footing, there- 
fore can spell out the spin-charge connections explicitly. 
From this spin-charge connection, we are able to classify 
all the possible excitations in the BLQH in a systematic 
way. The CB theory may also be applied to study the 
incoherent disordered insulating side and the quantum 
phase transitions between different ground states [13]. 
Using the CB approach, I also studied several interesting 
phenomena specific to im-balanced BLQH. Just like any 
Chern-Simon theory, the CB theory in BLQH has its own 
limitations: it is hard to incorporate the LLL projection 
( see however [14] ), some parameters can only be taken 
as phenomilogical parameters to be fitted into the micro- 
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scopic LLL+HF calculations or experimental data. It is 
an effective low energy theory, so special care is needed 
to capture some physics at microscopic length scales 

In this paper, we use both MCF approach and the CB 
approach developed in [12] to study spin polarized Tri- 
layer Quantum Hall (TLQH) systems at total filling fac- 
tor vt = 1. We also supplement the two approaches by 
wavefunction approach. TLQH is interesting from both 
experimental and theoretical sides. On the experimental 
side, TLQH systems have been fabricated in high mo- 
bility electron systems in the experimental group led by 
Shayegan [15]. On the theory side, because the three lay- 
ers play the roles of three flavors. The excitons in BLQH 
is the pairing of particle in one layer and the hole in the 
other. This pairing structure in BLQH is similar to the 
Cooper pairing in BCS superconductors up to a particle- 
hole transformation. While the excitons in TLQH is not 
obvious, because so far there is no analog of three flavors 
superconductors yet. It would be interesting to under- 
stand the structures of the excitons and all the possible 
excitations in TLQH. As shown in this paper, due to the 
particular structure pattern of the excitons, there are two 
coupled Goldstone modes in the TLQH ( see Eqn.33 ). In 
the presence of interlayer tunnelings ( with or without in- 
planc magnetic field ), the interference between the two 
gapless modes may also lead to many new phenomena 
which can not be seen in BLQH. 

There are at least two well defined regimes for TLQH. 
(I) Interlayer coherent regime: when the distance d bew- 
teen the two adjeent layers is sufficiently small, then all 
the three layers are strongly correlated. (II) Weakly- 
coupled regime: when d is sufficiently large, then all the 
three layers are are weakly coupled. Depending on the 
distance d, there could be many other possible regimes. 
In this paper, we focus on regime (I) which is the in- 
terlayer coherent regime. The system in the interlayer 
coherent regime (I) was discussed in [6,16] in LLL+ HF 
approach, two Goldstone modes are found. Here, we 
study the regime (I) in detail from both MCF and CB 
approaches and also stress the broken symmetry state in 
terms of ST/ (3) algebra. 

Consider a tri-layer system with N\, N2, N$ electrons 
in layer 1 ( the bottom layer ), layer 2 ( the middle layer 
) and layer 3 ( the top layer ) respectively in the presence 
of magnetic field B = V x A. We assume equal interlayer 
distance d between adjacent layers, the total filling fac- 
tor is vt — 1 and the spin is completely polarized. The 
Hamiltonian H = Hq + Hi nt is 

H = d'xcl(x) ^ c a (x) 

H lnt = ^ J d 2 xd 2 x'5p a {x)V a[j {x - x')5pp(x') (1) 

where electrons have bare mass m and carry charge — e; 
c a ,Sp a (x) = cJ t (x)c Q (x) — p a ,a = 1,2,3 are electron 



operators and normal ordered electron densities on the 
three layers. The intralayer interactions are Vn = V22 = 
V33 = e 2 /er, while interlayer interaction is V12 = V21 = 
V23 = V 32 = e 2 /eVr 2 + d 2 ,V 13 = V 31 = e 2 /eVr 2 +4d 2 
where e is the dielectric constant. 

The rest of the paper is organized as follows. In section 
II, we discuss MCF approach to TLQH and point out its 
weakness. In section III, we discuss CB approach and 
derive an effective action involving the charge and the 
two Goldstone modes. We also derive a dual action of 
the CB approach and make critical comparisons between 
this dual action and the MCF action. We classify all 
the possible excitations and discuss possible 2-body and 
3-body bound states and associated Kosterlize-Thouless 
transitions. In section IV, we will outline several salient 
features of correlated interlayer tunnelings with or with- 
out in-plane magnetic field in TLQH. We reach conclu- 
sions in the final section. At appropriate places in the 
paper, we point out the limitations of the CB approach. 

II. MUTUAL COMPOSITE FERMION 
APPROACH: 

We can extend the single-valued singular gauge trans- 
formation in [12] to tri-layer system: 

tj _ e iJ2 a pf d ^ X I d2x 'U a f)Pc,(S)arg(x~x')pp(x') ^ 

where the 3x3 symmetric matrix U and its inverse is 
given by: 

^=(101) U^ = \( 1-11] (3) 
\1 10/ z \ 1 11/ 

The Hamiltonian Eqn.l is transformed into: 

rr f,2 , f ,^ HfiV+fl(£)-fi5a(g)) 2 , _ . . 
Hq = Id X1pl{x) ^ ' ^ ' 

where the transformed fermion ip a (x) = UcafflU^ 1 is 
given by: 

lP a (x) = e i j d 2x 'U a ,[ar 9 (S-S')+arg(S'-S)]p,(3') c ^ ^ 

and the three Chern-Simon (CS) gauge fields a a in Eqn.4 
satisfies: V • a a = 0, V x a a = 2-KU a ppp(x) = 2tt[p{x) — 
p a {x)} where p(x) — J2 a Pa(%) 1S the total density of the 
system. 

In the following, we put Ti — c = e = e=l. At total 
filling factor vt = l,\7x A = 2irn where n = ni + n 2 + n 3 
is the total average electron density. By absorbing the av- 
erage values of CS gauge fields Vx < a a >— 2n(n — n a ) 
into the external gauge potential A* a = A— < a a >, we 
have: 
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H, 



(x) 



(-iV + A* a (x) - SS a (x)f 
2m 



1pa(x) (6) 



where V x A* a = 2irn a and V x 5S a = 2TtU a p8pp(x) — 
2n[Sp(x) — Sp a (x)] are the deviations from the corre- 
sponding average density. 

When d < d c /2, the strong intra- and inter-layer in- 
teractions renormalize the bare mass into two different 
effective masses m* = 7713,7712 [14]. MCF in each layer 
feel effective magnetic field £>* = V x A* a = 2Trn a , there- 
fore fill exactly one MCF Landau level. The energy gaps 
are simply the cyclotron gaps of the MCF Landau levels 

* -B a 

Transport properties: Adding three different source 

gauge potentials A 8 a to the three layers in Eqn.6, inte- 
grating out the MCF ip a first, then the mutual CS gauge 
fields, we can calculate the resistivity tensor by Kubo 
formula. We find that both Hall resistivity and two Hall 
drag resistivities are h/e 2 . This means if one drives cur- 
rent in one layer, the Hall resistivities in all the three 
layers are quantized at h/e 2 , while the longitudinal resis- 
tivities in all the three layers vanish. This is the hallmark 
of interlayer coherent quantum hall states. 

In the following, we confine to balanced case Ni = 
N 2 = N3 = N/3, imbalanced case can be discussed along 
the similar line developed in [12]. 

Fractional charges: Let's look at the charge of quasi- 
particles created by MCF field operators ^ a {x). If we in- 
sert one electron at one of the layers, say layer 1, from the 
singular gauge transformation in Eqn.2, we can see this 
is equivalent to insert one MCF in layer 1, at the same 
time, insert one flux quantum at layer 2 and another one 
at layer 3 directly above the position where we inserted 
the electron. The inserted flux quantum at layer 2 and 
layer 3 are in the opposite direction to the external mag- 
netic field, therefore promote one MCF in layer 2 and one 
in layer 3 to the second MCF Landau level in layer 2 and 
layer 3 respectively. Because inserting one electron leads 
to three MCFs in each layer, we conclude the charge of 
each MCF is -1/3. The above argument gives the correct 
total fractional charges of MCF, but it can not determine 
the relative charge distributions among the three layers. 
At mean field level, the energies of all the possible rela- 
tive charge differences are degenerate. The lowest energy 
configuration can only be determined by fluctuations. 

Left- and right- moving gapless modes: 5p a (x) can be 
expressed in terms of the CS gauge fields 5p a (x) = 

[U^^ap-^- where Uf 1 is given in Eqn.3. This constraint 
can be imposed by Lagrangian multipliers aP a which plays 
the role of time components of CS gauge fields. Integrat- 
ing out MCF tpi, -02, 1P3 to one-loop and carefully expand- 
ing the interlayer Coulomb interaction to the necessary 
order in the long-wavelength limit leads to an effective 
action for the three gauge fields. 
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9™' 9 [SeUj2 + {Sx+ ^ 
where a s = Si + 0,2 + S3, Si = S\ — S3, a r = Si + 03 — 2(?2 
which stand for center of mass, left and right channels re- 
spectively [17]. e c = (2ei + e 2 )/3,Xc = (2%i + X 2 )/3, e; = 
ei. Xl = Xi, £ r = (d + 2e 2 )/3, Xr = (xi + 2x 2 )/3, Se = 
£i — £2, fix — Xi ™ X2- The dielectric constants e a = 2 ™g, 
and the susceptibilities Xq = 2 mn* were calculated in 
single layer system in [18]. • • • are higher gradient terms 
and is the transverse component of gauge fields in 
Coulomb gauge V • S a = 0. The system is invariant 
under the Z 2 symmetry which exchanges layer 1 and 3, 
namely a c — > a c ,a; — > — a;,a r — > a r . This symmetry 
dictates that only a c a c , a c a r , a r a r , a; a; can appear in the 
Maxwell terms. The last two terms in Eqn. 7 which 
is the coupling between a c and a r can be shown to be 
irrelevant in the low energy limit. 

More intuitive choices are a u = a± — a 2 , a<j = a 2 — 03 
which is related to a/, a r by a; — a u + ad, a r = a u — ad- 
But a u and ad are coupled, the two normal modes are 
ai,a r instead of a u ,ad- As expected, there is a Chern- 
Simon term for the center of mass gauge field S s , while 
there are two Maxwell terms for the left and right chan- 
nel gauge fields Si, S r which are gapless. The two gapless 
modes stand for the relative charge density fluctuations 
among the three layers. After putting back K, c, e, e, we 
find the two spin-wave velocities in terms of experimen- 
tally measurable parameters: 



3(a; c *) 2 2ac d 
2irn 1 e A T- 



3«) 2 1 + 2 7 - 
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(8) 



2im l + 2 7 e 'H' v / 27ml + 27 

where we set oj* — lo* c1 and n is the total density, I is the 
magnetic length, 7 is the ratio of the two effective masses 
7 = m\/m\ and a ~ 1/137 is the fine structure constant. 

As in BLQH [12], we expect the second term domi- 
nates the first, then vi/v r <~ \fl + 27. namely, the ratio 
of the two velocities is determined by the ratio of the two 
effective masses. If we take the effective masses to be the 
band mass, then vi/v r <~ \/3- 

Topological excitations: As discussed in the previous 
paragraph, the MCF carries charge 1/3, but their relative 
charge distributions among the three layers are undeter- 
mined. They can be characterized by their (a c ,ai,a r ) 
charges (±1/3, qi, q r ) with qi,q r — 0, ±1, • • •. Exchang- 
ing a c leads to 1/r interaction between two MCF, while 
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exchanging ai , a r leads to logarithmic interactions which 
may lead to a 2-body bound state between two MCF with 
opposite (qi,q r )- The energy of this bound state with 
length L is A + + A_ + q 2 c ^ - ( qil q l2 + )hcu* InL/l 
where A + ,A_ are the core energies of QH and QP re- 
spectively. There are also possible 3-body bound states 
of three MCF with £. q u = £. 9™ = 0. 

Unfortunately, the gluing conditions ( or selection rules 
) of (q c ,qi,q r ) for realizable physical excitations arc not 
clear. This charge q c and two spin q r ) connections can 
only be easily established from CB approach presented 
in the next section. So we defer the discussion on the 
classification of all the possible excitations to the next 
section. 

Problems with MCF approach: 

All the criticisms on MCF approach to BLQH also hold 
to TLQH. In the following, we list just ones which are the 
most relevant to this paper. 

(a) It is easy to see that the spin wave dispersion in 
Eqn.8 remains linear u> ~ vk even in the d — > limit. 
This contradicts with the well established fact that in the 
d — > limit, the linear dispersion relation will be replaced 
by quadratic SU{3) Ferromagnetic spin- wave dispersion 
relation u) <~ k 2 due to the enlarged SU (3) symmetry at 
d^ 0. 

(b) The broken symmetry in the ground state is not 
obvious without resorting to the (111,111) wavefunction 
Eqn.28. The origin of the gapless mode is not clear. 

(c) The physical meaning of qi,q r is not clear. That 
they have to be integers was put in by hand in an ad hoc 
way. 

(d) The spin-charge connections in (q c ,qi,q r ) can not 
be extracted. 

(e) Even in the balanced case, there are two MCF cy- 
clotron gaps huj* a ,a = 1,2 at mean field theory. How- 
ever, there is only one charge gap in the system. It is not 
known how to reconcile this discrepancy within MCF ap- 
proach. 

(f) In the presence of interlayer tunneling, it is not 
known how to derive the tunneling term in a straightfor- 
ward way ( see section IV) . 

In the following, we will show that the alternative CB 
approach not only achieve all the results, but also can 
get rid of all these drawbacks. 

III. COMPOSITE BOSON APPROACH 

Composite boson approach was found to be much more 
effective than the MCF approach in BLQH. We expect 
it remains true in TLQH. In this section, we apply CB 
approach to study TLQH and find that it indeed avoids 
all the problems suffered in MCF approach listed in the 
last section. Instead of integrating out the charge de- 
gree of freedoms which was done in the previous CB ap- 
proach, we keep one charge sector and two spin sectors 



on the same footing and explicitly stress the spin and 
charge connections. We classify all the possible excita- 
tions, bound states and associated KT transitions. 

By defining the center of mass, left and right channels 
densities as: 

5p c = 5pi + Sp 2 + Sp 3 
Spi = Spi - Sp 3 

Sp r = Spi - 25p 2 + Sp3 (9) 
We can rewrite H int in Eqn.l as: 

Hint = ^SpiUSpx + ^Sp 2 USp 2 + ^Sp 3 USp 3 

+ 5piV5p 2 + Sp 2 V5p 3 + SpiWSps 

111 

= 2 S PcV c Sp c + -^SpiViSpt + -Sp r V r Sp r + Sp c V cr Sp r (10) 

where U = V n ,V = V 12 , W = V 13 and V c = f + ^ + 
- l(U-W),V r = I(^-^ + |), v cr = -Uv- 
W). The stability of the system dictates V c V r — V} r > 0. 
In the long wavelength limit qd <C 1, we only keep the 
leading terms: V c = V t = 2nd, V r = ^,V cr = ~^. 
Performing a singular gauge transformation: 

o (f) = e «/rfVar 3 (*-*')p(*') Ca(a?) (n) 

where p(x) = c\{x)c\{x) + cl(x)c 2 (x) + c' 3 (x)c 3 (x) is the 
total density of the tri-layer system. 

After absorbing the external gauge potential A into 
the Chern-Simon gauge potential a, we can transform the 
Hamiltonian Eqn.l into the Lagrangian in the Coulomb 
gauge: 

,t/^ ■ \, ,+ (— iTiV — ha(x)) 2 , 

C = (t>l(d T - lao)^ + 0t(f)i ^ 4>a{£) 

+ ia p+^ J d 2 x'6p c (x)V c (x - x')5p c (x') 

+ X - J d 2 x'5 Pl (x ) Vi (x - x')6 Pl {x') 

+ d 2 x'5p r (x)V r (x — x')5p r (x') 

+ / d 2 x' 8p c {x)V cr {x ~ x)5p r {x') - 7^-a (V x a) (12) 

J 27T 

In the Coulomb gauge, integrating out ao leads to the 
constraint: V x a = 2-n{^\(^ a — p). 

It can be shown that <j) a (x) satisfies all the boson com- 
mutation relations. We write the three bosons in terms 
of magnitude and phase 

<f> a - VPa + 5p a e t9 " (13) 

The boson commutation relations imply that 
[5p a (x), 4>b(x)] = ihS a bS(x — x'). For simplicity, in this 
paper, we only consider the balanced case, so we set 
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p a = p. The imbalanced case can be discussed along 
similar lines developed in [12]. 

We define the center of mass, left-moving and right- 
moving angles which are conjugate angle variables to the 
three densities defined in Eqn.9 



9c = @i + $2 + &3 

61 = 0i — #3 

9 r = 9i — 202 + 83 



(14) 



which satisfy the commutation relations [5p a (x), 6p{x' )] = 
A a ih6 a p5(x — x 1 ) where A = 3, 2, 6 for a = c,l, r. 

Substituting Eqnl3 into Eqn.12, neglecting the magni- 
tude fluctuations in the spatial gradient term which were 
shown to be irrelevant in BLQH in [12] and rewriting 
the spatial gradient term in terms of the three angles in 
Eqn.14, we find: 



c = is Pc (^d T e c -a ) + £ i [^ve c 



al 2 



l;Sp c V c (q)5p c 



-5p r d T 0T , 

3om 



(V9rY + -5p r Vr(q)Spr 



+ Sp c V cr (q)Sp r - — a (V x a) 
Ztt 



(15) 



In the balanced case, the symmetry is U(l) c x U(l)i x 
U(l) r x Z 2 where the first is a local gauge symmetry, 
the second and the third are global symmetries and the 
global Z2 symmetry is the exchange symmetry between 
layer 1 and layer 3. 

(a) Off-diagonal Algebraic order in the charge sector: 
At temperatures much lower than the vortex excitation 
energy, we can neglect vortex configurations in Eqn.15 
and only consider the low energy spin-wave excitation. 
The charge sector ( 8 C mode ) and the two spin sectors ( 
61, 8 r modes ) are essentially decoupled. It can be shown 
that the second to the last term 5p c V cr (q)5p r in Eqn.15 
which is the coupling between the charge sector and right- 
moving sector is irrelevant in the long wave-length limit. 
Therefore, the charge sector is essentially the same as the 
CSGL action in BLQH. Using the constraint a t — 27r ^ Pc , 
neglecting vortex excitations in the ground state and in- 
tegrating out 8p c leads to the effective action of 6 C : 

uj 2 +w% 



2 V 3 ; cV n Vc {q) + ^ 



hrWetfM (i6) 



where a; 2 - = u, 
frequency. 



§lQ 2 Vc(q) and u c = ^ is the cyclotron 
From Eqn.16, we can find the equal time correlator of 



— < 6 c (-q,-oj)6 c (q,uj) 



> 



= 3x^ + 0(1) 

q 2 q 



(17) 



which leads to the algebraic order: 



< e i(Oc(S)-O c m > = 



\x - y\ 3 



(18) 



We could define C = (6»i + 6 2 + 63) /3 = c /3, then the 
exponent will be 1/3. But when we consider topological 
vortex excitations, then e l9c may not be single valued. 
Therefore, 9 C is more fundamental than C . 

(b) Spin wave excitations: 

While in the spin sector, there are two neutral gapless 
modes: left-moving mode and right-moving mode [17]. 
Integrating out 5pi leads to 



2V l (q) y 2 



12m 



(19) 



where the dispersion relation of spin wave can be ex- 
tracted: 



^ = [^V l (q)]q 2 =vfq 2 



(20) 



In the long wave-length limit qd<g.l, the bare spin wave 
velocity is: 



2 P 47re j 
vf = — d 



m 3e 3me 
Integrating out Sp r leads to 



1 <ld T e r f + -^{^e r f 



2V r (q) y 6 



36m 



(21) 



(22) 



where the dispersion relation of spin wave can be ex- 
tracted: 



CO 



[— 

m 



V r (q)]q 



2 2 2 



(23) 



In the long wave-length limit qd <C 1, the bare spin wave 
velocity is: 



p Aire' 
m 9e 



d = v 2 /3 



(24) 



Both spin wave velocities should increase as the square 
root of the separation d when d < d c /2, their ratio 
vi/v r <~ \/3. At d = 0, vi = v r = 0. This is expected, 
because at d = the U(l)i x U(l) r x Z 2 symmetry is 
enlarged to SU(3)c, the spin wave of SU(3) isotropic 
ferromagnet uj ~ k 2 . 

We have also treated the coupling term between the 
charge and right-moving sector 6p c V cr (q)Sp r in Eqn.15 
exactly and found that in the q,u) — > limit, the mixed 
propagator takes the form: 



< e c (-q,-uj)6 r (q,u;) >^ 



2„2 



UJ < 



> 



(25) 
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which is dictated by < 9 r 9 r > instead of < 9 C 9 C >. Be- 
cause of the lu 2 prefactor, we can ignore < 9 c 9 r > in the 
low energy limit. 

We also found that Eqns.16, 17, 18 and Eqns.22, 23,24 
remain true in the long- wavelength and low energy limit, 
therefore confirm that Sp c V cr (q)Sp r is indeed irrelevant 
in the limit. 

As shown in BLQH in [12], the functional form of the 
spin sector in the CB theory is the same as EPQFM. 
Although there are no detailed microscopic calculations 
in TLQH yet except the preliminary work in [6,16], from 
the insights gained in BLQH, we can claim that the func- 
tional form of the spin sector in Eqn.15 ( or Eqn.33 ) 
achieved in the CB theory is correct. Again, it is hard to 
incorporate the LLL projection in this CB approach, So 
some parameters can not be evaluated within this ap- 
proach. For example, the two spin stiffnesses for the 
left and right-moving sectors in Eqn.15 ( or or Eqn.33 
) should be completely determined by the Coulomb in- 
teractions instead of being dependent of the band mass, 
they may be estimated by the microscopic LLL+ HF ap- 
proach. However, the ratio of the two velocities in Eqn.24 
is independent of the band mass, it is interesting to see if 
this ratio \/3 remains correct in the microscopic LLL+HF 
estimation. Just like in BLQH, we can simply take p s i, Vi 
and p sr ,V r in Eqn.15 as phenomilogical parameters to 
be fitted into the LLL + HF estimations. However, just 
like in BLQH [10,12,19], because the qualitative correct 
ground state wavefunction in TLQH is still unknown, it 
is still difficult to estimate these parameters with con- 
trolled accuracy even in the LLL+HF approach. We can 
simply fit these parameters into experimental data. 

(c) Dual action: 

Just like in BLQH discussed in [12], we can per- 
form a duality transformation on Eqn.15 to obtain a 
dual action in terms of three vortex currents J?'' c ' l ' r = 
^e fJ ,uxdudx9 c ,Lr and three dual gauge fields b^b^b 7 ^: 

2ir 

C d = -inbfavxdvbl - iA^e^xd^ + i—b^J^ 

777 1 -* — » 

+ - d ob c a ) 2 + x W(9)(V x b c ) 

- iA^e^xdvbi + mb 1 ^ 1 
+ ^(dX - d,b l a f + 1(V x ?)VH«7)(V x ?) 
-iA r s ^xd v b r x + i^j; r 

+ ^(d a b r - d b r a ) 2 + l(v x b r )v r (q)(v x v) 

+ (V x b c )V cr (q)(W x b r ) (26) 

where the three source gauge fields are A c s = A\ stl + 
A2.su A^sui A s ^ = A\ s ^ — A^s^, A s ^ = Axgjj, — 2A2.su + 

We can compare this dual action derived from CB ap- 



proach with Eqn.7 derived in MCF approach. Just like 
in BLQH, the three topological vortex currents are miss- 
ing in Eqn.7, while Xc XhXr, Sx terms are extra spurious 
terms which break SU (3) symmetry even in the d — > 
limit. If we drop all the x terms, use bare mass and 
add the three vortex currents terms in Eqn.7, the result- 
ing effective action will be identical to the dual action 
Eqn.26. Neglecting the vortex currents, we can identify 
the two spin wave velocities easily from the dual action 
which are the same as found previously from Eqn.15. It 
is also easy to show that the last term in Eqn.26 which is 
the coupling between the charge and right-moving sector 
is irrelevant in the low energy limit. We conclude that 
Eqns.15 and its dual action Eqn. 26 are the correct and 
complete effective actions. 
(d) Topological excitations: 

Any topological excitations are characterized by three 
winding numbers A9i = 27rmi,A0 2 = 27rm 2 ,A0 3 = 
27rm3, or equivalently, A9 C — 2ir(mi + + m^) = 
2irm c ,A9i = 27t(toi— to 3 ) = 2nmi,A9 r = 27r(mi — 2to2 + 
to 3 ) = 2nm r . It is important to stress that the three fun- 
damental angles are #1, #2, 6> 3 instead of 9 c ,9i,9 r . There- 
fore, mi,TO2,m 3 are three independent integers, while 
to c , mi,m r are not. 

There are following 6 kinds of fundamental topological 
excitations: A<f>i — ±2n or A0 2 = ±27r or A</> 3 — ±2ir, 
namely (mi,ra 2 ,m 3 ) = (±1, 0, 0), (0, ±1, 0), (0, 0, ±1). 
They correspond to inserting one flux quantum in layer 1 
or 2 or 3 in the same or opposite direction as the external 
magnetic field. Let's classify the topological excitations 
in terms of (g, m;, m r ) where q is the fractional charge of 
the topological excitations in the following table. 





(±1,0,0) 


(0,±1,0) 


(0,0, ±1) 


m c 


±1 


±1 


±1 


mi 


±1 





±1 


m r 


±1 


±2 


±1 


q 


±1/3 


±1/3 


±1/3 



Table 1: The fractional charge in balanced TLQH 

The fractional charges in the table were determined 
from the constraint V x a = 2ir5p c and the finiteness of 
the energy in the charge sector: 

q= -^-(j)a-dl = ^- x ^j>\76 c -dl = ^m c (27) 

There are the following two possible bound states: 
(1) Three Charge neutral two-body bound states: 

• Layer 1: (±1/3, ±1, ±1) with energy E x = E c+ + 
E c- - m + 2lr (P sl + P°r) In ^ when R > R c > I. 

• Layer 2: (±1/3, 0, +2) with energy E 2 = E c+ + 
E c _--^ + 2ir{4p sr )\n^-. 

• Layer 3: (±1/3, ±1,±1) with energy E$ — E\ 
which is dictated by Z 2 symmetry. 



6 



The corresponding 2-body Kosterlize-Thouless (KT) 
transitions are k B T^ T1 = hsT^^ = %(p s i + p sr ) which 
is dictated by Z 2 symmetry and k B T"^ T2 = 2-Kp sr above 
which the QP and QH pairs are liberated into free ones. 
If setting p s i — 3 p sr , then there is only one 2-body KT 
temperature k B T^ T = 2irp sr = As explained previ- 
ously, microscopic calculations in the LLL [22] are needed 
to estimate the two spin stiffnesses p s i and p sr . In gen- 
eral, there should be two 2-body KT transition temper- 
atures k B Tf T1 = k B T^ b T3 ^ k B Tf T2 . The spin wave 
excitations will turn into these Charge neutral two-body 
bound states at large wavevector ( or short distance ). 
The two-body bound state behaves as a boson. 
(1) Two Charge ±1 three-body bound states: 
The three-body bound states consists of three quasi- 
particles (±1/3, ±1, ±1), (±1/3, 0, t2), (±1/3, ±1) 
located at the three corners of a triangle ( Fig.l). 



implies mi = 771,2 = WI3 = ffi and 77i c = 3m. Eqn.27 im- 
plies the charge q — \m c = m must be an integer. This 
proof rigorously rules out the possibility of the existence 
of deconfined fractional charges. We conclude that any 
deconfined charge must have an integral charge, m = 1 
corresponds to inserting one flux quantum through all 
the three layers which is conventional charge 1 excita- 
tion. Splitting the whole flux quantum into three fluxes 
penetrating the three layers at three different positions 
will turn into the three-body bound state with the same 
charge shown in Fig.l. It is still not known if the energy 
of this conventional charge 1 excitation is lower than the 
three body bound state. 



IV. BROKEN SYMMETRY 



(1/3,1,1) 




(1/3,0,-2) R 23 (1/3,-1,1) 

Fig 1: The lowest energy charged excitation is a three-body 
bound state of three ±1/3 charged quasi-particles sitting on the 
three corners of a triangle with R12 = R23 dictated by the Z2 sym- 
metry. The three quantum numbers (q,mi,m r ) are enclosed in the 
parenthesis 

Its energy E 3b = 3E C± + ^(-L_ + _L_ + _L_) + 
2tt [2p sr In ^ + (p s i - p sr ) In ^ + 2p sr In ^f] when R > 

R c 3> I. Minimizing E 3 b with respect to R12, R13, R23 

2 

leads to two optimal separations: R \2 — R023 — — 

OD7T Psr 

which is dictated by the Z2 symmetry and i? i3 = 
18tt(p^ _ p r . If setting p s i = 3p sr , then the three quasi- 
particles are located on the corners of a isosceles with 

2 

length R = 36 ^ p , the corresponding 3-body Kosterlize- 
Thouless (KT) transition is the same as the two-body 
bound state k B TfP T = k B Tf T = 2-Kp sr = f^. These 
charged excitations behave as fermions and are the main 
sources of dissipations in transport experiments. In gen- 
eral, Tj^ T should be different from the two 2-body KT 
transition temperatures. Again, microscopic calculations 
of the two spin stiffnesses p s i and p sr in the LLL [22] are 
needed to determine the three KT transition tempera- 
tures. 

Let's look at the interesting possibility of deconfined 
( or free ) 1/2 charged excitations. Because any excita- 
tions with non-vanishing mi or m r will be confined, so 
any deconfined excitations must have mi = m r = which 



The broken symmetry and associated Goldstone modes 
can be easily seen from wavefunction approach. In the 
first quantization, in the d — > limit, the ground state 
trial wavefunction is the (111, 111) state: 

JVi N 2 JVi N 3 N 2 N 3 

*ni,in = nn^ ~ "jjini^ ~ " ii ii ' " • 

i—lj—l 

JVi N 2 N 3 

X Y[( Ul - Uj) Y[(Vi - Vj) Y[(Wi - Wj) (28) 

i<j i<j i<j 

where u,v,w are the coordinates in layer 1, 2 and 3 re- 
spectively. 

In the second quantization, a trial wavefunction was 
written in [16]: 

M-l 

i*> = n 7! (^<i+c I t „, 2 +^ci, 3 )io> 

M-l 

= IT ^={l + e W -Cl A C m .2 + e^Cl^C m .2) 

M-l 

x I] <2|0> (29) 

m=0 

where 6 U = 6>i - 2 ,0 d = 63 - 9 2 and M = N is the 
angular momentum quantum number corresponding to 
the edge. We can interpret 6 U mode in Eqn.29 as a "up" 
pairing between an electron in layer 1 and a hole in layer 
2 which leads to a 9 U molecule, similarly, Od mode as a " 
down " pairing between an electron in layer 3 and a hole 
in layer 2 which leads to a 9d molecule ( Fig. 2). 
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is easier to perform transport experiments in the counter- 
now channel than in the right-moving channel. 




(3a) 



Fig 2: Two pairing modes of TLQH. Solid dots are electrons, 
open dots are holes. 

Note that when 9 U = 9d = 0, I* > is the ground 
state \G >, it can be shown [19] that its projection on 
(Ni,N 2 , N 3 ) sector is (111, 111) wavefunction in Eqn.28. 
Note that Eqn.28 is correct only in d — > limit where 
the symmetry is enlarged to SU(3). In the SU(3) sym- 
metric case, the pairing between any two layers of the 
three layers is equally likely. However, at any finite d, 
Eqn.28 and Eqn.29 may even not be qualitatively cor- 
rect [19], because the SU(3) symmetry is reduced to 
U(l) x U(l) x U(l) x Z 2 symmetry ( where the three 
17(1) symmetry correspond to c a (x) — > e l0a c a (x), a 
1 . 2, 3, Ci <-» c 3 ) The pairings between any two layers out 
of the three layers are not equivalent any more. The two 
lowest energy pairing modes are shown in Fig. 2. 

SU (3) has 8 generators which can be written in terms 
of c a : 

Si(x) = c\(x)d(x) - cl(x)c 3 (x) 

S r (x) = c\(x)ci(x) + cl(x)c 3 (x) - 2c\(x)c 2 (x) 
S\ 2 (x) = c\(x)c 2 (x), S 2 i{x) = S\ 2 (x) = c\(x)ci(x) 
S 23 (x) = c\(x)c 3 (x), S 32 (x) = Sl 3 (x) = c\{x)c 2 (x) 




21- 



(3b) 



Fig 3: (3a) Left-moving superfluid channel which is the counter- 
flow channel (3b) Right-moving superfluid channel [17] 



V. CORRELATED INTERLAYER TUNNELINGS 

So far, we assume that the tunnclings are tuned to 
be zero. In this section, we will outline several salient 
features in the tunnclings of TLQH. We assume the tun- 
neling amplitude from layer 1 to layer 2 is the same as 
that of from layer 2 to layer 3. The tunneling term is: 



H t = tc\c 2 + h.c. + tc\c 3 + h.c. — t(j)\(j} 2 + h.c. + t(j)\cj) 3 + h.c. 

(31) 



6 

tp(cos 9 U + cos 6d) = 2tp cos — cos 



Si 3 (x) = c[(x)c 3 (x), S 3 i{x) 



S\ 3 (x) 



where Si , S r are the two ( left and right ) generators in 
Cartan subalgebra ( Note that the algebra is not SU(2) in 
spin j = 1 representation which is also 3 dimensional as 
used in [6] ). In the balanced case N\ — N 2 = N 3 = N/3, 
we have < Si >=< S r >= 0. The ground state in Eqn.28 
breaks two of the three U(l) symmetries, therefore there 
are two Goldstone modes in the corresponding two order 
parameters < S i2 >= fe l6 ^,< S 23 >= fe"* 9 ". 

When 6 U molecule and 0d molecule in Fig. 2 move in 
the opposite direction, the currents in layer 2 cancel each 
other, it leads to the left-moving superfluid channel 0i ( 
Fig. 3a ). While when 8 U molecule and 9d molecule move 
in the same direction, the currents in layer 2 add and flow 
in the opposite direction to the currents in layer 1 and 
layer 3, it leads to the right-moving superfluid channel 
9 r ( Fig. 3b ). In BLQH, there is only one superfluid 
channel which is the counter- flow channel [5]. In TLQH, 
there are two superfluid channels which are left-moving 
and right-moving channels, the left-moving channel cor- 
responds to the counter-flow channel [17]. Obviously, it 



where 6 U = 6± - 6 2 ,6 d = 6 3 - 6 2 
c\{x)c\{x) (30) Integrating out the gapped charge sector 6 C leads to 
the effective tunneling action: 



c = llxi(drOi) 2 + Ps^Qif] + l -[ Xr {d T e r ) 2 + Psr {ve r f] 



+ 2ipcos — cos — 
' 2 2 



(32) 



where, as noted previously, the two spin susceptibilities 
XhXr and the two spin-stiffnesses p s i 7 p S r are phenomi- 
logical parameters to be fitted into the LLL+HF calcu- 
lations or experimental data. 

Unfortunately, the above action is not very useful 
for any practical calculations, because the two angles 
0i = 6 u -6 d = 0i-0sA = 0u + 0d = 9 1 -2d 2 + d 3 are not 
independent angles when considering the vortex excita- 
tions. We can see this fact by noting that 0i + 9 r = 29 u 
or 81 — 9 r = —29d has to be twice an angle which is a 
constraint between 9i and 9 U . While 9 U and 9d are two 
independent angles. So it is much more useful to write 
the above tunneling action in terms of the two indepen- 
dent angles 6 u ,6a: 
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C=^( X l+ Xr)(d T O u ) 2 + \{psl + Psr){V6 u f + tp COS 6 U 
+ \{Xl + Xr){d T O d f + ^(p sl + p sr ){V6 d ) 2 + tp COS 6 d 

+ (xi - xr)d T u d T 9 d + ( Psl - Psr )ve u ■ ve d (33) 

where there is a Z 2 symmetry under 9 U <-» 9 d . The 
tunneling currents in three layers can be expressed in 
terms of the two independent angles: I\ = tpsm6 u , I 2 = 
—tp(sm6 u + sm6d),h = tpsm8 d . They satisfy the cur- 
rent conservation I\ + I 2 + I3 = 0. 

From Eqn.33, we can see there is a coupling between 
U and 9 d . If there were not this coupling, we would have 
two independent BLQH tunneling actions. It is this cou- 
pling which lead to the correlated tunnclings of TLQH. 
For example, the current out of layer 1 Ii = tp sin 6 U only 
depends on 6 U explicitly, while 6 U couples to 6 d . It was 
shown that the single gapless mode in BLQH leads to a 
sharp zero-bias peak in the interlayer tunneling of BLQH, 
the experimentally observed dissipations can only come 
from external sources such as disorder [21]. If the cou- 
pling in TLQH even in the absence of external sources 
will lead to a broad zero-bias tunneling peak in TLQH 
will be investigated in [22]. 

In the presence of in-plane magnetic field B\\ — 
(B X7 B y ), the gauge invariance dictates the tunneling 
term: 

H t = tp(cos(6 u -Q-x) + cos(9 d + Q-x)) (34) 



In BLQH, it was found that the in-plane field will split 
the zero-bias peak, the splitting gives a direct measure- 
ment of the spin-wave velocity of the single Goldstonc 
mode in BLQH [21]. However, in TLQH, there are two 
coupled Goldstonc modes with two different velocities 
vi,v r . In [22], we will study how the zero-bias peak shifts 
and how to relate the shift to the two spin wave velocities 
or their ratio in the presence of in-plane magnetic field. 

We can also see that if there were not the coupling 
between 6 U and 9 d , we would have two independent 
Pokrovsky-Talapov (PT) models. In BLQH which is de- 
scribed by a single PT model, it was found that when 
the applied in-plane magnetic field is larger than a criti- 
cal field B > By , there is a phase transition from a com- 
mensurate state to an incommensurate state (C-IC) with 
broken translational symmetry. When B > BT,, there 
is a finite temperature KT transition which restores the 
translation symmetry by means of dislocations in the do- 
main wall structure in the incommensurate phase [9,11]. 
In TLQH which is described by the above two coupled 
PT model [20] , several kinds of C-IC transitions and as- 
sociated several kinds of KT transitions are expected and 
will be studied in a separate publication [22]. 



VI. CONCLUSIONS: 



In this paper, we study the interlayer coherent incom- 
pressible phase in spin polarized Trilayer Quantum Hall 
systems at total filling factor vt = 1 from three ap- 
proaches: MCF, CB and wavefunction approach. Just 
like in TLQH, CB approach is superior than MCF ap- 
proach. The Hall and Hall drag resistivities are found 
to be quantized at h/e 2 . Two neutral gapless modes 
which are left and right moving modes with linear dis- 
persion relations vi/ r = LOi/ r k are identified and the ratio 
of the two velocities is close to \/3. The novel excita- 
tion spectra are classified into two classes: (1) Charge 
neutral bosonic two-body bound states. (3) Charge ±1 
fermionic three-body bound states. In general, there are 
two 2-body KT transition temperatures and one 3-body 
KT transition temperature. Microscopic calculations in 
the LLL [22] arc needed to roughly estimate the three 
KT transition temperatures. The Charge ±1 three-body 
bound states shown in Fig. 1 may be the main dissipation 
source of the interlayer tunneling. The broken symme- 
try in terms of SU{2>) algebra is studied. The effective 
action of the two coupled Goldstone modes are derived. 
The excitonic structure in TLQH was shown in Fig. 2. 
In BLQH, there is only one superfluid channel which is 
the counter-flow channel [5]. In TLQH, there are two su- 
perfluid channels which are left-moving and right-moving 
channels ( Fig. 3 ), the left-moving channel corresponds 
to the counter-flow channel ( Fig. 3a ). Obviously, it is 
easier to perform transport experiments in the counter- 
flow channel than in the right-moving channel. In the 
presence of interlayer tunnclings, the coupling between 
the two Goldstone modes may lead to the broadening 
of the zero-bias tunneling peak, in contrast to the zero- 
bias peak in BLQH. The precise nature of the broadening 
and several other possible new phenomena unique to the 
TLQH in the presence of in-plane magnetic field will be 
studied in a forthcoming publication. 

We expect that the effective action of the two coupled 
Goldstone modes and the qualitative physical pictures 
of TLQH achieved from CB approach in this paper are 
correct. However, the parameters in the effective action 
need to be estimated from microscopic calculations in the 
LLL [22] . It is interesting to see if the ratio vi jv r <~ \/3 
still holds in the microscopic calculations. As stressed 
in [10,19] for BLQH, the trial wavefunction Eqn.28 and 
Eqn.29 may not even be qualitatively correct in TLQH, 
so it is still difficult to estimate the ratio with controlled 
accuracy even in the LLL+HF approach. Eventually, the 
ratio need to be measured by experiments. 

I thank J. K. Jain for helpful discussions. 
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